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Abstract 

For a two-dimensional quantum mechanical problem, we obtain a generalized power- 
series expansion of the S'-matrix that can be done near an arbitrary point on the Riemann 
surface of the energy, similarly to the standard effective range expansion. In order to do 
this, we consider the Jost-function and analytically factorize its momentum dependence 
that causes the Jost function to be a multi-valued function. The remaining single-valued 
function of the energy is then expanded in the power-series near an arbitrary point in 
the complex energy plane. A systematic and accurate procedure has been developed for 
calculating the expansion coefficients. This makes it possible to obtain a semi-analytic 
expression for the Jost-function (and therefore for the S-matrix) near an arbitrary point 
on the Riemann surface and use it, for example, to locate the spectral points (bound 
and resonant states) as the S-matrix poles. The method is applied to a model simlar to 
those used in the theory of quantum dots. 
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1 Introduction 



It is probable that an average reader of this journal perceives the one- and two-dimensional 
problems as simplified toy models of quantum mechanics. Although such an attitude has its 
roots in the standard courses of quantum mechanics, this perception is far from being ade- 
quate. First of all, modern nano-technology allows us to fabricate the microscopic quantum 
devices that behave and can be described as one- or two-dimensional. The two-dimensional 
tunneling of particles play an important role in superconductive tunnel junctions and even in 
some biological molecules [Ij. Besides that, the corresponding quantum mechanical problems 
are not mathematically simple as one may think. Indeed, in contrast to the motion of a particle 
in three-dimensional space, the one-dimensional motion of the same particle on an infinite line 
is inherently a two-channel problem, where the channels are the left and right halves of the 
line (see, for example, Refs. [2,3J). As far as the two-dimensional scattering is concerned, its 
amplitude as a function of the energy has not a square-root but a logarithmic branching point 
(see Refs. [4]-[l0] as well as the subsequent sections of the present paper). Therefore from 
both pure mathematical and practical points of view, the two-dimensional quantum problem 
is worthwhile to consider. 

In the present paper, we focus on solving the two-dimensional problem with the help of power 
series that are similar to famous effective- range expansion (where hk and 5o are the collision 
momentum and the S-wave phase-shift) 



in terms of the so called scattering length a, effective radius tq, etc., introduced long ago in 
nuclear physics [TT]. 

As we already mentioned, the energy dependent functions of the two-dimensional problem have 
a logarithmic branching point at the threshold. As a result there is a controversy concerning 
the two-dimensional analog of Eq. ([1]). Some authors [6]-[8l[T2] define the scattering length a' 
by including it in the logarithmic term. 



(here 7 is the Euler's constant) while the others [4l[5j try to preserve the traditional form of 
the right hand side of Eq. ([1]) and move the logarithmic term to the left hand side. 



k cot So{k) 
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where In a' = —ixjo!'. 

We look at this problem from a more general point of view. What is actually done in the 
original effective range expansion ([1]) is the constructing of the function R{E) = kcot5o{k) 
in which the square-root branching of A; ~ ±^/E at the threshold is compensated by exactly 
the same branching of 5o{k) ~ ±\/E. As a result the function R{E) does not have branching 
points and is a single-valued analytic function of the energy E ~ k'^, and therefore can be 
expanded in a convergent series R{E) = + o-iE + 02-^^ + • ■ ■ , which is given by Eq. ([T]). 

From this reasoning a next logical step immediately follows: the function R{E) can be ex- 
panded in a more general power series R{E) = bQ + bi(E — Eq) + b2{E — Eq)"^ + ■ ■ ■ around an 
arbitrary complex energy Eq within the domain of its analyticity. In Refs. [T21II1], we realized 
this idea for the three-dimensional single-channel and multi-channel problems. In doing this, 
instead of using R{E), we expanded the analytic single-valued parts of the Jost functions 
(or Jost matrices in the multi-channel case) after explicit separation of the factors that are 
responsible for the branching. 

In the present paper, we use the same approach as in Refs. [T3l[T4] to obtain similar expansions 
of the Jost functions for the two-dimensional problem. First, we analyze the analytic struc- 
ture of the Jost functions and split them in the single-valued and logarithmically branching 
parts. Then, we derive a set of differential equations that determine the single-valued parts. 
And finally, we look for the solutions of these equations in the form of power series of the 
energy. The series and ([3]) together with simple recipes for calculating any number of their 
expansion coefficients, are easily obtained from our more general expansions that are done 
around an arbitrary complex point Eq. Using two-dimensional model potential related 
to quantum dot theory, we numerically demonstrate the efficiency and accuracy of the 
proposed method. 

2 Jost function 

Radial part ui of the wave function describing the motion of a particle of mass jj, with the 
energy E \n dt circularly symmetric two-dimensional potential U{r) = h^V{r)/2ii obeys the 
differential equation (a review of the partial-wave analysis for the two-dimensional scattering 
is given in the Appendix [Al) 

«.,(£;,r)=0, (4) 



ar r 
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where ^ = A + 1/2 is the angular momentum and A = —1/2, 1/2, 3/2, .... At large distances, 
where V{r) — i- 0, the radial equation simplifies 



A(A + 1) 



ue{E, r) ^ , r — )■ oo . (5) 



This is the Riccati-Bessel equation. As its two linearly independent solutions, we can choose 
either the Riccat-Bessel and Riccati-Neumann functions jx{kr) and y\{kr), or the two Riccati- 
Hankel functions h^^\kr). Any other solution of Eq. (E]) is a superposition of the two linearly 
independent solutions. In particular, we can write the asymptotics of the physical wave function 
as a linear combination of the Riccati-Hankel functions, 

u,{E,r) ft\E)ht\^,{kr) + ft'''\E)h[%{kr) , (6) 

where the energy-dependent combination coefficients f^"^°^^\E) are called the Jost functions. 
When r — )■ oo, the Riccati-Hankel functions represent the incoming and outgoing circular 
waves. Indeed, 



' V 2 |2|— !>00 

where H^(^\z) are the cylindrical Hankel functions. The Jost functions f\™^°^^\E) are there- 
fore the asymptotic amplitudes of the incoming and outgoing waves. Since the flux of the 
particles is conserving, for real E we have \ff'^\E)\ = \f}°^^\E)\. Actually, these two func- 
tions are related to each other at different complex values of E as well. Some of such symmetry 
properties can be established using the semi-analytic structure of them that is derived in the 
subsequent sections. It can also be shown that the partial wave S'-matrix is the ratio of these 
functions 

ft\E) 

and that zeros of the Jost function f^™\E) are the discrete spectral points 8n, 



^^(E) = , (8) 



ft\^n) = 0, (9) 
i.e. the bound and resonant states of the system. 
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3 Transformation of the radial equation 



Our goal is to establish the analytic structure of the Jost function, i.e. to find such an ex- 
pression for it where all possible nonanalytic dependencies on the energy are given explicitly. 
This can be done if we transform the second-order radial equation (jlj) into an equivalent set 
of first-order equations. 



The transformation is done using a method which is known in the theory of differential equa- 
tions as the variation parameters method [TSldFj- Following this method, we look for the 
unknown function ui{E,r) in the form similar to its asymptotics (E]), but with the combina- 
tion coefficients being new unknown functions of r, 

u,{E,r) = Ft\E,r)h[z\/,{kr) + Ft'''\E,r)h[%{kr) , (10) 

where Ff^/°"*)(E , r) are the new unknown functions. Since instead of one unknown function, 
we introduce two of them, they cannot be independent of each other. In principle, we can 
impose any reasonable condition relating them. Looking at the asymptotics (El), we see that 

/f H^) = lim Ff HE,r) , ft'''\E) = lim Ff"*^(E,r) . (11) 

r— ^-oo r— >oo 

Therefore at large distances our new functions become constants, and we should have 



drF, 



(in) 



drF, 



(out) 



/iW(A;r) = . 



(12) 



As the additional condition imposed on these functions, we demand that the relation (IT^ is 
valid not only at large r but at all distances. In the variation parameters method this con- 
dition is known as the Lagrange condition. In our case this condition makes F^^^^°^^\E,r) 
to be the Jost functions for the potential which is cut-off at the radius r (in the spirit of the 
variable-phase approach). 



Substituting the ansatz (1101) into the radial equation ([4]) and using the Lagrange condition 
(11^ together with known Wronskian of the Riccati-Hankel functions, 

h[z\kr)drh^^\kr) - h'^^\kr)drh{z\kr) = 2ik , (13) 

we obtain a set of two first-order equations which are equivalent to the original radial equation 



drF, 



(in) 



drF, 



(out) 



2ik 



2ik ^ 



(14) 
(15) 



The boundary conditions for these equations follow from the requirement that the wave func- 
tion must be regular everywhere. In particular, this means that U£{E,0) = 0. It seems that 
this is not the case because both h^^\kr) and h{^\kr) that are present in the expression 
( ITOl) . are singular at r = 0. Their singularities, however, can cancel each other. 



(16) 



if they are superimposed with a same coefficient. This can be achieved if both Fj^^^\E,r) and 



F^°^^\E,r) have the same value at r = 0, 



n(out) 



{E,0) 



Their common value at r = determines the overall normalization of the wave function and 
therefore can be chosen arbitrarily. To be consistent, we chose it to be 1/2, which makes 
U£{E,r) to behave near the origin exactly as the Riccati-Bessel function and thus the solution 
with the boundary conditions 



is what is called the regular solution in the theory of three-dimensional scattering. 



(out) , 



1 



(17) 



For our goal of expressing the non-analytic dependencies of the Jost functions in an explicit 
form, it is more convenient to re-write the ansatz (fTUf) in terms of the Riccati-Bessel and 
Riccati-Neumann functions, 

ut{E, r) = A,{E, r)jx{kr) - B,{E, r)yx{kr) , (18) 

and to obtain the corresponding equations for the unknown functions Ae{E,r) and B^^E^r). 
Since 

h^^\z)=3y^{z)±iyy^{z) , (19) 
this is most simply achieved by making the following linear combinations of Eqs. (IT4l [T5|) 

- Ft\E,r) + Ft'''\E,r) , 



ME,r) 
B,{E,r) 



This gives 



drBi 



Fr\E,r)-Fr\E, 



-^yxV{Aijx - Beyx) 



-p\V{Aijx - Biyx) 



(20) 
(21) 

(22) 
(23) 
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with the boundary conditions 

A,(E,0) = 1, 5,(E,0)=0. (24) 

4 Complex rotation 

Suppose that the potential V{r) is cut off at certain radius r = R, then the right-hand sides 
of the sets of equations (11411150 or (12211230 vanish at r > _R and thus the derivatives on the 
left-hand sides of these equations become zero, i.e. the functions or and Bi do 

not change beyond this point. Therefore, in the spirit of the variable phase approach, the 
functions F^^^^"^^\E,r) are the Jost functions for the potential which is cut off at the point r. 
Generally speaking, when the potential asymptotically vanishes at large distances, we should 
expect the convergence of the limits (ITTj) . 

Therefore, the Jost functions can be calculated by numerical integration of the differential 
equations (1 1411 151) or (I22II23I) from r = up to a sufficiently large radius R where the limits 
(ITT!) are reached within a required accuracy. This works perfectly for real values of the energy 
E. However, when we consider complex energies (for example, in search for resonances), 
a technical difficulty arises. This difficulty is caused by the asymptotic behavior ([7]) of the 
Riccati-Hankel functions. 

When k is complex, either h^^\kr) or h^^\kr) exponentially diverges, depending on the sign 
of ImA;. As a result, either the first or the second of the equations (11411150 does not give a 
numerically convergent solution. This difficulty is circumvented by using the deformed integra- 
tion path shown in figure [H Instead of integrating the differential equations along the real axis 
from r = to r = i?, we can reach the final point via the intermediate point r = i?' in the 
complex plane. Moreover, we can safely ignore the arc R'R since the potential is practically 
zero at that distance. 

Such a complex rotation helps because the asymptotic behavior (divergent or convergent) of 
the functions hf'\kr) is determined by the sign of ImA;. For a given k = \k\e^'^, we can 
always find such a rotation angle in r = |r|e~*^ that the product kr = \kr\e^^'^~^^ has either 
positive or negative (or even zero) imaginary part. Various technical details of using complex 
rotation in calculating the Jost functions and Jost matrices can be found in P1I31 IT31II41ITBI - I25] . 
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5 Explicit separation of the non-analytic factors 

In order to establish the analytic structure of the Jost functions, we need to have a closer look at 
the structure of the Riccati-Bessel and Riccat-Neumann functions. The following expressions 
for them (they can be derived using formulae 9.1.2, 9.1.10, and 9.1.11 of Ref. [IT]) are the 
most useful for this 



where 



n=0 

oo 

y^ikr) = k-'Y,k''^gi'\r) + h{k)j^{k 

n=0 



n!r(n + A + 3/2) V2 
If A is integer then the expansion of y\{kr) is also simple 

^" ^ ' n!r(n-A + l/2) V2 
h{k) = . 

However, for a half-integer A, we have a more difficult case: 

{\-n - 1/2)! /r\2"-A 



for any A . 



(25) 
(26) 

(27) 

(28) 
(29) 



A/vrn! 



^ n ^ A 



2 ' 



(30) 



■1) 



n X 5 |^/;(^ 



+ 1) + ipin - X + I)] /rN 2n-A 

2. 



v^n!(n- A - i)! 



A + I ^ < oo 



h{k) = ^ln(^) , (31) 

where R is an arbitrary number (in the units of length). It is arbitrary because any increase or 
decrease of h{k) cased by the change of i? is compensated by the corresponding change in the 
first term of Eq. ( !30l) . The parameter R is introduced to separate the r and k dependencies 
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in the original term containing \n{kr/2) = ln(/ci?/2) + \n{r/R) and to have dimensionless 
products under the logarithm. In practical calculations the parameter R can always be taken 
as the unit of the length, i.e. R= 1. The ^/^-function in Eq. (!30f) is defined as follows [17j 



ijj{n) 




-7 



7 + E 



m=l 



n-1 



m 



-1 



n ^ 2 , 



n = 1 , 



(32) 



where 7 = 0.577. . . is the Euler constant. 

Eqs. (!25l [261) represent the Riccati-Bessel and Riccati-Neuman functions in the form of infi- 
nite series. Each term of these series is a product of a function depending on k and another 
function depending on r, i.e. the k and r dependencies are given in a separable form. 

What do the above formulae tell us about the Jost functions? The functions j\{kr) and y\{kr) 
are involved in the coefficients of the differential equations (12211230 that determine the Jost 
functions. This means that the Jost functions are not single valued functions of the energy. 
Indeed, for each choice of E, we have two possible values of the momentum 



The index A = —1/2, 1/2, 3/2, . . . of the Riccati functions is a half-integer. This means that 
the differential equations involve such multi-valued functions as square-root and logarithm of 
the momentum. 

Therefore the Jost functions are defined on a complicated Riemann surface and the threshold 
point = is a branching point of this surface. It would be desirable to find an expression 
for the Jost functions in terms of the powers of \/k, the logarithmic function h{k), and some 
entire single valued functions of E. In order to do this, we notice that the series in Eqs. 
(I25II26I) involve only even powers of k, i.e. the powers of the energy k'^ = 2jj,E/h^. Since for 
any finite r these series are absolutely and uniformly convergent on the whole complex plane 
of E, they define some entire functions, i.e. 




yx{kr) 



k-'UE,r) + k^+'h{k)jx{E,r) , 



(33) 



(34) 
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where the "tilded" functions 

UE,r) = f:(^-^ye\r), (35) 

n=0 ^ ^ 

UE,r) = f:(^)%i^Hr), (36) 

n=0 ^ ^ 

are single-valued entire functions of complex variable E. 

Let us find a similar structure for the functions Ai{E,r) and Bi{E,r) and through them for 
^^(m/out)^^^ r). For this, we replace the set of equations (12211230 with their linear combinations. 
Namely, we multiply Eq. ( !23l) by h(k) and subtract the result from Eq. ( !22l ): and as the second 
equation, we take Eq. ( !23l) multiplied by k^^'^'^^^\ As a result, we obtain: 

dr{Ae-hB,) = -\{yx-hjy)V{A,j^-B,yy) , (37) 

= -A;-2(^+i)j,\/(A,j, - B,y,) . (38) 
Now, taking into account Eqs. ([331 l34j), we see that 

yx - hjx = k-^yx + k^^^hjx - k^+^hjx = k'^y^ 

and 

Aejx-Beyx = A^k^+'jx - Bek'^'y^ - Bek^+'hj^ 

= k^+\Ae~hBe)jx-k-^Beyx. 
Substituting these expressions into Eqs. (1371 [38|) . we have 



dr{Ae~hB,) = -k-'^'+''>yxV[k^+\Ae-hBe)jx~k-^B,yx] , (39) 

drk-^''+''>Be = -k-^'+''>jxV[k'+\Ae-hBe)jx-k-'B,yx\ . (40) 
If we introduce the "tilded" functions 

AeiE,r) = A,{E ,r) ~ h{k)B,{E ,k) , (41) 

B,{E,r) = k-^^^+^~^B,{E,r) , (42) 
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then Eqs. (13^ |4D|) assume the following form 

drAe = -yxV [Aijx-Beyx) , (43) 

drBe = -jxV(^Aax-Beyx) (44) 
with the boundary conditions [that follow from i^^] 

ME,0) = 1, Be{E,0)=0. (45) 

For any finite r, all the coefficient functions in Eqs. (1431 |44I) are entire functions of the 
parameter E and the boundary conditions are i?-independent. According to the Poincare 
theorem \76] the solutions of these equations, i.e. the functions Ai(E,r) and Bi(E,r), are 
entire (analytic single-valued) functions of the complex variable E. 

Therefore the structure we wanted to find is as follows: 

Ae{E,r) = A,{E,r) + k^^+^h{k)B,{E,r) , (46) 

B,{E,r) = k^^+^Bt{E,r) , (47) 

where Ai{E,r) and Bi{E,r) are single-valued analytic functions of E. Apart from these 
single-valued functions, the original functions Ag and Be involve the factors k^^^^ and h{k). 
Since A is half-integer, the power (2A + 1) is always even and thus k"^^^^ is also single-valued 
function of E, but h{k) has a logarithmic branching point at = 0. 

The functions i^^*^™/""*) have similar structure 

Ft\E,r) = ^{Ae-^Be) = ^{ME,r) + e'+'[h{k)-^]Be{E,r)] , (48) 

Ft'''\E,r) = ^{Ae + tBe) = ^{ME,r) + k''+'[h{k)+t]Be{E,r)y (49) 

6 Analytic structure of the Jost functions 

What we have established in the previous Section, is valid for any complex E and any finite 
distance r. In other words, so far we have established that the Jost functions (ITIi) have the 
structure (1451 |4^ if the potential is cut off at certain radius r = R (does not matter how 
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large R is). The problem is that to prove analyticity of Ai{E,r) and Bi{E^r) with respect to 
variable E, we used the Poincare theorem which requires that the coefficients of Eqs. (I43ll44p 
be holomorphic functions of E. If is a real positive number and the potential is of a short 
range, this is true even if r — )■ oo. Indeed, in such a case both jx{kr) and y\{kr) oscillate 
with finite amplitudes even at infinity and thus the coefficients of Eqs. ( I43B41) simply tend 
to zero (i.e. remain holomorphic) when r — )■ oo. If however E is negative or complex, then 
generally speaking this is not true. As we will see shortly, there still is a domain of complex E 
where the coefficients remain holomorphic. In other words, if we extend r to infinity, we have 
to narrow the domain of E. 



The Riccati-Bessel and Riccati-Neumann functions are linear combinations of the Riccat- 
Hankel functions and thus at large distances behave as exponential functions ([7]). If the mo- 
mentum has a nonzero imaginary part, then one or the other of these exponentials is diverging 
and thus both j\{kr) and y\{kr) tend to infinity when r — )• oo. To some extent the situation 
can be saved by using a short-range (exponentially decaying) potential V{r) ~ exp(— r^r), 
which compensates the divergence of j\{kr) and y\{kr) within certain domain V of the 
complex E'-plane along its real axis. The borders of the domain V are determined by the 
requirement that none of the coefficients of Eqs. (14311440 are divergent. The behavior (con- 
vergent or divergent) of these coefficients is determined by the product exp(±2iA;r) exp(— ?7r). 
For a given rj it is not difficult to find the domain V, 



V= {E: 

which gives the condition 



2Im ^2iiE/n? <r]^ , (50) 



(Im Ef <-^ + ^Re E . (51) 



^2 ^ fi^ri^ h^'rf 
16/12 + ^ 

Similar analyticity domain was obtained by Motovilov [27j for the three-dimensional multi 
channel T-matrix, using a rigorous analysis of the corresponding scattering operators. 



The faster the potential decays, the wider is the domain. An example of such a domain is 
shown in figure [2] for the model used in section It is a parabolic domain along the real 
axis, whose border is shown by the solid curve. It crosses the real axis at ~ —0.0332 (in 
donor Hartree units). When r — )• oo, we can only use the Poincare theorem within V, and 
thus we can say that at least within this domain the functions ^^(-E', oo) and Bi{E, oo) are 
the holomorphic functions of variable E. 

If the potential V{r) (or at least its long-range tail) is an analytic function of complex variable 
r and exponentially decays along any ray r = |r|e*^ within certain sector of the complex r- 
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plane, then the domain T> can be extended by using the complex rotation described in section 
m With a complex radius, the product exp(±2iA;r) exp(— ?7r) vanishes at infinity if E is within 

V = {E : |2Im(A;r)| < r/Rer} , (52) 

which generalizes Eq. ( ISOll . It is easy to show that \^ E = |ii^|exp(ix) then such a domain 
can be defined by the following inequality 

sm — h < ; — ; — . 53 

V2 / 8/i|E| ^ ' 

With ^ = this condition is transformed into (ISTl) . An example of such a domain with the 
rotation angle 6 = 0.05n for the potential (!84l) is shown by the dashed curve in Fig. [21 

The physically interesting domain of the i^-plane where the structure (1451 1491) . can be used in 
practical calculations, lies on the positive real axis (scattering) and in the close vicinity below 
it (pronounced resonances). Therefore, we can say that to all practical purposes this structure 
is valid at an arbitrary point E. 

If we denote the asymptotic values (within V) of the "tilded" functions as 

de{E) = lim A,{E,r) , k{E) = lim 5,(E,r) , (54) 

r-->oo r—^ca 

Then the Jost functions and the S-matrix can be written as follows 

/f)(E) = ^{~aeiE) + k''^'[hik)-^UE)] , (55) 

= l{ME) + k''^'[h{k)+z]k{E)] , (56) 

^ ^ ~ae{E) + k'^-^'[h{k)+tME) 

a,{E) + k^^+^[h{k) - i]k{E) ■ 

To find the Jost functions or the 5'-matrix on any sheet of the Riemann surface, we need 
to calculate the functions ai{E) and b£{E) only once (because they are single valued). The 
choice of the sheet is determined by an appropriate choice of the value of the logarithmic 
function h{k). Please note that k"^^^^ = k'^^ is a single-valued function of E. 
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7 Power-series expansions of the Jost functions 



The functions ai(E) and b^^E) are holomorphic (i.e. single-valued and analytic) and therefore 
can be expanded in Taylor series near any point Eq within the domain V of the complex energy 
plane. The expansion around the point Eq = will give us the standard effective- range series. 
But we can also do such an expansion near an arbitrary point, 

oo 

de{E) = J2^n\Eo){E - EoT , (58) 

n=0 

oo 

k{E) = J2l^^n\Eo){E-Eor . (59) 



n=0 



How the expansion coefficients an^ and f3n^ can be found? For this purpose, we can derive 
differential equations, the solutions of which asymptotically tend to «„ and /3n ■ Indeed, 
such an expansion can be done at any fixed radius r because the functions Ai(E,r) and 
Bi{E,r) reach their limits (!54f) at r if the potential is cut off at this radius (in the spirit of 
the variable-phase approach). Therefore for each r, we have 

oo 

A, {E,t) = Y,^'n\Eo,r){E-Eor , (60) 

n=0 

oo 

B, iE,r) = Y,Bl'\Eo,r){E-Eor , (61) 



n=0 



where 



a 



(Eo) = lim ^(f)(Eo, r) , /3i'\Eo) = lim B^^\Eo, t) . (62) 



Therefore, the differential equations mentioned above, should determine the functions ^^^(£"0, r) 
and Bn\EQ,r). In order to obtain such equations, we expand the "tilded" functions jA(-E', ?^) 
and yx{E,r) in the Taylor series near an arbitary point Eq 

00 

UE,r) = ^.W(Eo,r)(E-Eo)" , (63) 

n=0 
00 

UE,r) = J2ci'\Eo,r){E~Eor , (64) 



ra=0 
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which are more general expansions than the series (13511361) for the particular case of the thresh- 
old energy Eq = 0. Any number of the expansion coefficients si^'\Eo,r) and ch\Eo,r) can 
be found using the recurrence relations derived in the Appendix [B] 

Substituting the expansions ( [5D1I61II63II541) into Eqs. (14311440 . and equalizing the factors of the 
same powers of (E — Eq), we obtain the equations we are looking for, 

SrAi? = - E cl^'V (4%™ - 8f><^>) , (65) 

i+j+k=n 

= - E ■'S"'^ i^P^' - 4'') . (66) 

i+j+k=n 

with the boundary conditions 

^f(Eo,0) = 5„o, i3f(^o,0) = 0, n = 0,1,2,3,... (67) 

These conditions follow from the fact that the corresponding boundary conditions (1450 do not 
depend on E. Therefore, starting with the initial values (1571) at r = 0, and numerically solving 
first A^ + 1 pairs of differential equations of the system (!55l up to a sufficiently large radius 
'"max, we obtain first + 1 expansion coefficients 

a^^\Eo) = (^0, rmax) , /3f (^o) = i^f (^o, r„.ax) , n = 0, 1, 2, . . . , iV (68) 
These coefficients give us the following approximate formulae for the Jost functions 

1 ^ 

ft\E) ^ 2 S + k^'+Mk) - i]P^\E,)] {E - EoT , (69) 

n=0 
1 ^ 

ft"\E) ^ 2 5Z i^niEo) + e'^'m + ^H'\Eo)} {E - E^T , (70) 

n=0 

which are valid for any complex value of E within a domain around the chosen central point 
Eq. Apparently, the closer E is to Eq, the better is the accuracy of these formulae. It is 
interesting to note that Eq. (9) of Ref. ^lOj is the first term of our Eq. ( 15^ for the particular 
case of i?o = 0. 

An alternative way of using formulae ( I5WU11 is to treat the expansion coefficients an\EQ), 
(3n\Eo), n = 0,1,..., N as fitting parameters. Adjusting them in such a way that the 
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corresponding cross section (see Appendix IA.4p reproduces experimental data in vicinity of 
a real energy E^, one then can use the Jost function (15^ at the nearby complex energies 
for locating possible resonances. The obvious advantage of such an approach is that the 
resonance energy and the width are deduced directly from experimental data using correct 
analytic structure of the S'-matrix. 

8 Effective- range expansion 

Far away from the interaction region the radial wave function (1181) is a linear combination of 
the Riccati-Bessel and Riccati-Neumann functions 

ueiE,r) ^ aeiE)jxikr)-be{E)yx{kr) , (71) 

where 

ae{E) = lim ME,r) , 6,(E) = lim Be{E,r) . (72) 

r— >oo r— >oo 

The functions jx and yx in (iTii) can be written in their asymptotic form, 

jx{kr) — > sin ( kr - ^ ) , (73) 

r— >-oo \ Z J 

yx{kr) — > -cosf/cr-^j , (74) 



1 — ^oo 



which gives 



ue{E, r) — )■ ai{E) sin ( kr — + be{E) cos (kr — ^ 



A^sin 



kr-^ + 6,{E) 



where ag and bg are replaced with their common normalization factor and the scattering 
phase shift 5e, 

ae(E) = N cos 5i{E) , (75) 

bi{E) = N sin 6e{E) . (76) 
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Using the relations (jUl |47]) at large distances (r — )■ oo), 

a,{E) = ~a,{E) + e'+'h{k)k{E) , (77) 

be{E) = e'+%{E), (78) 

we can construct the so called effective-range function which is a holomorphic function of the 
energy. This is done by taking the ratio 

cot Of 



be k^^+^be 

and moving all the "troublesome" terms and factors which may generate singularities, to the 
left hand side of the equation, 

be 

[cot 6eiE) - h{k)] = [cot 6e{E) - h{k)\ = . (79) 

beiE) 

Both the numerator and denominator in the last ratio can be written in the form of power 
series §^ M) with = 

oo 

k''+' [cot SeiE) - hik)] = ^ . (80) 

Y.f^n\0)E- 

n=0 

Using Eq. (3.6.22) of the book by Abramowitz et al., 

ao + aix + a2x'^ -{ Oq T f ai bA f a2 6i(ai6o - ao&i) 2 

bQ + bix + b2x'^ ^ 60 L \aQ bo J \ao aob^ bo J 

the division of two polynomials in Eq. ( !80l) is done as follows 

-1 (e) 

[cot Se{E) - h{k)] = -^ + + ■ ■ ■ , (81) 
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where the scattering length a*^^^ and the effective radius Tq^-* for the state with the angular 
momentum ^ are given by 

On 



9 A numerical example related to quantum dot the- 
ory 

To demonstarate how the proposed method works, we use the following circularly-symmetric 
potential, which is motivated by the models that are currently used in the theory of quantum 
dots, 

U{r) = Voir -r 0)6-'/^ , (84) 

with Vq = 25, tq = 2, and R = 2, where Vq (as well as all the energies in this example) is 
measured in the so called "donor Hartree units" and the distances in the units of "donor Bohr 
radius" , which were chosen to be 10.96 meV and 101.89 A, respectively. These values for the 
units are relevant to the motion of electrons in the semiconductor material GaAs [28], where 
the effective electron mass \s fi = 0.063me (with mg being free electron mass). 

Although, strictly, the potential ( !84l) should be considered as an abstract quantum-mechanical 
"toy" model, we chose its shape in such a way that it resembles the potentials that are cur- 
rently used to describe two-dimensional quantum dots (see, for example, Refs. [28-31J). As 
is seen in FigjH our potential has a repulsive barrier which is not present in the traditional 
quantum-dot models. The main reason for introducing such a barrier was to enrich our "toy" 
model spectrum with resonances. However, one can argue that such a barrier may appear in 
real quantum dots as well. Indeed, when electrons fill up the lower levels of a dot, they should 
repel each other and tend to stay mostly at its periphery. This means that for an additional 
incoming electron the attractive force at the centre is reduced and a repulsion appears at the 
border. In other words, the original empty-dot confining-potential (shown with the dashed 
curve) is transformed into something that looks like our "toy" potential. Of course, this spec- 
ulative reasoning does not mean that we claim that our potential is anything more than an 
abstract model. 
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Since nothing special is associated with the angular momentum, we only consider here the 
S-wave states {i = 0). For such a case, the potential (!54l) supports three bound states and 
a squence of resonances. These spectral points (given in Tableland shown in Fig. [3]) were 
located using the exact approach, i.e. as the roots of Eq. ([9]), where /j"'' is the asymptotic 
value dm) of the solution of Eq. (fMl) . 

To make sure that we did not miss any of the bound states and/or narrow resonances, we 
calculated the S'-wave scattering phase-shift and checked if it obeys the Levinson's theorem. In 
Refs. ^[32l[33] it was shown that in the absence of a zero-energy bound state for the P-wave 
and always for the S'-wave, this theorem is the same as for the three-dimensional scattering, 
namely, 

5e{0) - 5e{oo) = ttN, , (85) 

where A^^ is the number of bound states with the angular momentum i. If the energy moves 
to the right along the real axis, the phase shift increases by vr near each resonance which is 
not far from the real axis. The smaller is the width, the more sharp is the increase. When 
calculating the phase shift numerically, it is easy to miss a sharp jump corresponding a narrow 
resonance. The curve "A" in Fig. [5] is an example of such omissions (the first two resonances 
are missed because of a too large step along the E'-axis). The correct phase shift is shown 
by the curve " B" . It starts with Stt at the threshold and tends to zero at the infinity, in 
accordance with Eq. (!55l) . 

Calculating the first two expansion coefficients and using Eqs. (I82ll83p . we found the following 
scattering length and effective radius, 

ao = -0.4521260323 [dimensionless] , tq = 0.0586790752 [length^] . 

As a first test of the expansions (I69ll70p . we performed them at several scattering energies 
(i.e. on the real energy axis) and compared the approximate cross section obtained from the 
approximate Jost functions (see Appendix IA.4I) with the corresponding exact cross section 
that was calculated using the exact Jost functions via numerical integration of the system of 
differential equations (11411150 . Fig. [5] shows the exact cross section in the interval E E (0, 10] 
(thick curve) and the approximate cross sections (thin curves) when only the first five terms 
of the series (I69II70I) were taken into account for Eq being 1, 5, and 7. It is seen that within 
rather wide interval around each Eq the expansion reproduces the cross section very well even 
with all its zigzags. 

The next step was to test our expansions at complex energies. To begin with, we performed 
them around a point on the real axis, namely, around Eq = 7 (far away from the threshold 
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energy) and looked at the Jost function at the nearby complex energies. Why 7? Simply 
because there is a resonance not far from this point (third resonance of Table [1]). To check 
the accuracy of the expansion, we compared the approximate values of f^^^\E) at various 
points around Eq with the corresponding exact values of the Jost function. Apparently, the 
closer the point E is to the center of the expansion, the more accurate should be the result. 
Fig. [7] shows three closed contours around £"0 = 7. Within the smallest of them the relative 
error of f^™\E) obtained by the expansion (15^ with = 4 is less than 1%. The other 
two contours show the domains of 5% and 10% accuracy. The important fact is that even if 
the expansion is done on the real axis, the semi-analytic formulae (I69ll70p remain valid at the 
nearby complex points. 

The star in Fig. [7]is a resonant zero of the exact Jost function. As is seen, the 1%-contour has a 
"dent" near this point. The reason for it is that in calculating the relative error, we have an ex- 
act value of f^^^\E) in the denominator and this value is zero at the resonance. By the way, the 
approximate Jost function with = 4 has zero at E = 7.1051679246 - |0. 5683685515 
which is not far from its exact position. This means that the expansion done on the real axis 
can be used for locating narrow resonances. 

Finally, we tested the expansion around a point in the fourth quadrant (where the resonances 
are) of the complex energy plane. When solving the differential equations (I65II65P we used the 
complex rotation of the coordinate (see Secjlj) with such an angle 9 that Im (kor) = (where 
ko is the momentum corresponding to Eq). This guarantees that Eq is within the domain V 
(see SecE]). 

Fig. |8] shows the exact positions of two resonances (indicated with stars), the center of the 
expansion (cross) at Eq = 7.55 — 21.06 which is in the middle between them, and two pairs of 
the approximate locations of these resonances: open circles for three terms of the expansion 
and filled circles for five expansion terms. It is seen that the expansion converges, i.e. the 
more terms are taken into account, the more accurately the resonances are reproduced. It 
should be noted that the chosen position of Eq is the "worst case" . If we move Eq a bit closer 
towards one of the resonances, it is reproduced much more accurately. 



10 Conclusion 

in this paper, we show that the Jost function for the two-dimension scattering can be written 
as a sum of two terms, one of which is an analytic single-valued function of the energy E 
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while the other term can be factorized in an analytic function of E and a logarithmic func- 
tion of the momentum. This means that the (logarithmic) branching point of the Riemann 
energy-surface is given in the Jost function explicitly via the logarithmic factor. The remain- 
ing energy-dependent functions are defined on single energy plane which does not have any 
branching points anymore. For these energy-dependent functions, we derive a system of first- 
order differential equations. Then, using the fact that the functions are analytic within certain 
domain V, we expand them in the power series around an arbitrary point Eq E V and obtain 
a system of differential equations that determine the expansion coefficients. 

A systematic procedure developed in this paper, allows us to accurately calculate the power 
series expansion of the Jost function practically at any point on the Riemann surface of the 
energy. Actually, the expansion is done for the single-valued functions of the energy, while the 
choice of the sheet of the Riemann surface is done by appropriately choosing the sheet of the 
logarithmic function of the momentum. 

The method suggested in this paper, makes it possible to obtain a semi-analytic expression 
for the two-dimensional Jost function (and therefore for the corresponding ^-matrix) near an 
arbitrary point on the Riemann surface and thus to locate the resonant states as the S-matrix 
poles. Alternatively, the expansion can be used to parametrize experimental data, where the 
unknown expansion coefficients are the fitting parameters. Such a parametrization will have the 
correct analytic structure. After fitting the data given at real energies, one can use the semi- 
analytic Jost function to search for resonances in the nearby domain of the Riemann surface. 
The efficiency and accuracy of the suggested expansion is demonstrated by an example of a 
two-dimensional model potential. 

APPENDICES 

A Two-dimensional partial-wave decomposition 

The partial-wave decomposition of the wave function, scattering amplitude, and cross section 
for a particle moving on a plane, is done using the cylindrical coordinates where the z-ax\s 
(perpendicular to the plane) is needed to define the orbital angular momentum. All the 
steps of such a decomposition are similar to the three-dimensional case, but the resulting 
formulae are not obvious and cannot be easily obtained from the corresponding 3-D-analysis. 
The derivations of various formulae of this type are given in several different papers (see, for 
example, Refs. [MIES]). Usually these derivations are very concise with many details omitted. 
Since such derivations are not present in the standard textbooks on quantum mechanics, we 
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feel that it is worthwhile to collect everything in one place. This is why we include this 
Appendix. 



A.l Radial Schrodinger equation 

Consider a particle of mass moving on a plane and being affected by a force that is described 
by a potential U{r), which is assumed to be of a short-range and circularly symmetric, 

U{f) = U{\f\) . 

To make the derivations simple, we assume that our particle does not have spin (this restriction 
can be easily revoked later). In the coordinate representation, the Hamiltonian H of such a 
particle is most conveniently expressed using the polar coordinates, 

r dr \ dr ) h^r^ ' 



where 



(9^ 



is the two-dimensional operator of the square of the angular momentum. Its eigenfunctions 
3^m(<^) obeying the equation 



and normalized as 



are easy to find. 



[ yMym'{^)dip = 6mm' , (86) 

yM^^e''^'^, m = 0,±l,±2,... (87) 

V ZTT 



From the definition of the Fourier series on the interval [0, 27r] it follows that 

+00 
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Each value of the angular momentum (except for zero) is represented twice: with two oppo- 
site signs. Classically, these two states correspond to the motion of the particle at the same 
distance r from the center and with the same velocity, but at different sides of the center (see 

Fig. m- 

The quantum number of the angular momentum i = \m\ is always non-negative and irrespec- 
tive of its magnitude (if i ^ 0) the vector i can have two (only two) directions: up or down 
(like the spin 1/2). The quantum number m = ±i\s its 2;-component. In principle, we can use 
the same notation for the eigenfunctions of the operator as in the three-dimensional case, 
namely, with two subscripts. However, because of the relation m = ±i, the subscript i 
is redundant. 



This can be formulated in a different way. The functions (1870 form a complete ortho-normal 
set on the interval ip E [0,27r]. This means that any (reasonable) function /((/?) defined on 
this interval, can be written as their linear combination, and such a combination can be written 
in the following two (equivalent) ways 



+0O oo 

m=—OQ £=0 m=±i 

In other words the following summations are equivalent 



^|m|m 



E --EE 



m=— oo £=0 m=±£ 

The operator corresponding to the quantum number m is obtained as follows. The gradient 
operator in the cylindrical coordinates is 

V = r— + ip-— + z— , (89 
or r Oifi oz 

where 

f = [0 X z\ , = [z X r\ , z = [f X if] 
are the corresponding unit vectors. Then 

-> ^ hr f 2,1 d ^ (9 

[r X pi = rr X —V = — z-— ip— 

I I \ r dip Oz 

and thus 

h d 



■'Z 



i dip 
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Apparently, both and iz commute with the Hamiltonian. The quantum numbers £ and 
m are therefore conserving. When specifying m, we implicitly specify the quantum number £ 
as well. This means that the quantum state of the particle is determined by two conserving 
quantum numbers, namely, the energy E and the 2;-component of the angular momentum m. 
The corresponding wave function, obeying the Schrodinger equation. 

Hi) Emir) = EipEmir) , 

can be factorized in the radial and angular parts 



VEm[r) = J= ymi'f) 



(90) 



where y/r in the denominator is introduced to obtain the radial equation without the first 
derivative. Substituting this factorized form into the Schrodinger equation, we obtain 



m 



1/4 



- V{r) 



Um{E,r) 







where k is the wave number (linear momentum) defined by 

2/i 



E 



(91) 



(92) 



and V{r) is the reduced (in the units of [length] ^) potential 

V{r) = ^J^U{r) . 

Noting that Eq. (19 ID is exactly the same for both choices of the sign for m, we conclude that 
UmiE,r) actually depends on i but not on m. The radial equation can therefore be re-written 
in the way we used to see it in the three-dimensional problems 



Ue{E,r) 







where we introduced 



X = i- 



(93) 



(94) 



and did the replacement 



1 
4 



+ 



A(A + 1) 
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Formally, Eq. (193 looks exactly like the radial equation of the three-dimensional problem. 
The only difference is that A is not an integer number 

_1 1 3 5 
~ ^2' 2' 2' 2' 

This simple fact makes a huge difference: it changes the analytic properties of the Jost function 
and thus the S'-matrix, because the Riccat-Neumann function y\{kr) with a half-integer A has 
a logarithmic branching point on the Riemann surface of the energy |T7] . 



A. 2 Plane- wave and circular waves 

Consider a two-dimensional plane wave normalized to the (5-function 



m = ^ , {k'\k) = 5{k'~k)= \5{k' - k)5{v' - ^) , (95) 

where ^9 is the polar angle of the momentum hk. This plane wave can be expanded over the 
full set {y} of the angular functions (!57l) . 

Y,cim{kr)y^{^) , (96) 



gjfcr gjfcr cos ^p 



27r 27r 

im 



where the x-axis is directed along the coordinate vector r. The expansion coefficients 



can be found using the integral representation of the Bessel function [IT] 

27ri™ /n 



e' 



2n 



iiz cos ^-mv) ^ _ ^i{-z cos ^-m^) _ ^gg^ 



2712'" 27r Jo 

Comparing Eq. ( !97l) with ( !98il and using the symmetry property of the Bessel function 

Jm{—z) = (— l)™Jm(-2), we See that 

am{kr) = —=Jm{kr) (99) 
V 27r 
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and thus 

ikr -I 1 1 +°° 

Using another symmetry property, J^m{z) = {—l)'^Jm{z), we see that the product i"^Jm{kr) 
does not depend on the sign of m and thus this expansion can be re-written as 

oo 



27r 



1 

2^ 



^ oo 

— e^z^ cos{£(f) Ji{kr) 



(101) 



where is the "multiplicity" of an £-state, i.e. is the analog of the factor {2£ + 1) of the 
3D-case, 

' 1 , i = 0, 

2 , i>0 . 



(102) 



Expressing the Bessel function via the Riccati-Hankel functions, 



JAz) 



2'KZ 



1-1/2^ 



we obtain the following decomposition of the 2D plane wave in the incoming (— ) and outgoing 
(+) circular waves 



27r 27rvfcr 



where A is defined by Eq. ((54 



(103) 



=0 ■m=±l 



In the above, we assumed that vector r was directed along the x-axis. If this is not the case, 
then the dot-product 

kr = kr{cos ipk cos ipr + sin ipk sin ipr) = kr cos{ipr — fk) 

depends on the two polar angles. In this general case the plane wave is expanded over two sets 
of functions {3^(v5fe)} and {y{^r)} depending on the angles of the momentum and co-ordinate 
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vectors. In a similar way as we did it above, it is not difficult to show that 



e 



-CX3 1=0 m=±e 

^ oo 

— ^ e^z^ cos [i{Lpr - ipk)] Ji{kr) 

-i= \h[-\kr) + h^^\kr)] y^^Myu^k) 

(E,r)3^;(y,,)3^„(^,) , (104) 



where 



uf\E,r) = h^~\kr) + h^^\kr) = 2jx{kr) (105) 

is a regular solution of the radial Schrodinger equation ( !93l ) for the case V{r) = 0. These 
partial-wave decompositions can be conveniently written in the following symbolic form 

1^) = \kim)ym{^k) , \^=Y1 I^M3^m(</^r) , (106) 

£m £m 



{rim\ki'm') = 6u>6mm'i'^Mkr) = 5u'5mm'i\l ^^ix{kr) , (107) 



{f\kim) = -l-j,{kr)y*^i^r) , {k\rim) = (-^Y ^ -l-j,{kr)y*^icp,) , (108) 
m\k'e'm') = l6{k - k')6fj'6^m' , {rim\r'i'm') = -6{r - r')6a'6mm' , (109) 

K T 

Y\k^m){k^m\kdk = l , / ^ |r£m)(r£m|r rfr = 1 . (110) 
A. 3 Scattering wave function 

The plane wave (11040 is a scattering wave function '(/'g(r) for the particular case of V{r) = 0. 
Apparently, the structure of its partial-wave decomposition should be the same for all potentials 

= ^=Y.''u,{E,r)y*^{^r)ym{^k) , 
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where the factor is determined by the choice of the potential and the collision energy (for 
the free motion, = 1 at all energies). The purpose of this factor is to always have exactly 
the same normalization, namely. 



m^k') = s[k-k') . (Ill) 

An appropriate value for can be found as follows. The Riccati-Hankel functions h^^\kr) are 
two linearly independent solutions of the radial Schrodinger equation (!^ without the potential 
term. This means that for a short-range potential, its solution asymptotically behaves as a 
linear combination of the Riccati-Hankel functions, 

Ue{E,r) ft'\E)h[-\kr) + f^'''''\E)h[^\kr) , (112) 

where the combination coefficients depend on the energy and are called the Jost functions. In 
fact, they are the amplitudes of the incoming (— ) and outgoing (+) circular waves. On the 
other hand, at large distances the wave function V^g(r) consists of the two parts: the initial 
(incident) wave 'ip^^\'r) and a scattered circular wave that goes in all directions with certain 
amplitude F, 



pikr 



If ?/'^(r) is properly normalized then ip^^\r) = e*^''/(27r). This means that the radial wave 
function at large distances should also be split in two parts one of which coincides with the 
function (IIOSD . In doing such a splitting of the function (I112p . we obtain 



u,{E,r) /f) 

r— j-oo 

and see that = l/ft\ i-e- 

V2^rfr\E) 



(out) 

h 



^.^r-) = ___^-^^^z%(E,r)3^;(^,)3^m(^.) ■ (113) 



A. 4 Cross section 

Defining the partial-wave S'-matrix and the amplitude, 

28 



and using 

as well as the fact that 



-iexp [i{kr — X7c/2)] 





1 


1 


2^ 




1 




2^ 


2n 


cos 



COS [i{(pk - (fr)] 



= 
> 



we can write the asymptotic behaviour of the scattering wave function as 



r— !>oo Zn 



Akr 



(114) 



where ip = ipk — (pr is the scattering angle and the total scattering amplitude has the following 
partial-wave expansion 



(115) 



1=0 



The cross section for a two-dimensional scattering has the units of length. The number of 
particles scattered into the angle spanned by the arc r dip, is the product of the flux in that 
radial direction and the length of the arc. The corresponding cross section da is defined as 
such a length that after its multiplication by the total incoming flux, it gives the same number 
of particles, i.e. 



•(in) 



da = ji°'''\ip)rdip 



Using standard definition for the particle flux j = h/ (2i/i)(^*V?/' — with the operator 

V given by Eq. (15^ . it is not difficult to find that the incoming [corresponding to the first 
term of the wave function (11141) ] and outgoing (obtained from the second term of the same 
wave function) fluxes are 



3 



■(in) 



nk 



■ (out) 



(27r)2;x ' 

and thus the differential cross section is 



hk\j^{E,^)[ 

(27r)V 
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Using the integral 



/ cosiiip) cos{e' cp)dip ^ I 27r , ^ = f = > 



27r 



the total cross section can be written as follows 



/.27r 

cr ^ \T{E,ip)fdip = y2'^i^ 

Jo 



a, = 27re,|/,(£;)|' = I -1|' , 



where cr^ is the partial-wave cross section. 



(116) 



B Expansion coefficients for the holomorphic parts of 
the Riccati-Bessel and Riccati-Neumann functions 

The Riccati-Bessel and Riccati-Neumann functions j\{kr) and y\{kr) can be written in the 
following factorized form 



yx{kr) = k-%{E,r) + k^^'h{k)jx{E,r) , 



A+l, 



(117) 

(118) 



where the "tilded" functions are holomorphic with respect to the energy variable E. This 
means that we can expand them in the Taylor series, 



UE,r) = J2'n\Eo,r)iE-Eoy 



n=0 



n=0 



near an arbitrary point Eq. The expansion coefficients, 

si'\Eo,r) = ^^h{E,r) 



ci'\Eo,r) 



n\ dE"" 

1 d"" 
rJ.dE^ 

30 



y\{E,r] 



E=Eo 



E=Eq 



(119) 
(120) 

(121) 

(122) 



are expressed via the corresponding derivatives. In order to find them, we notice that 

h^k"^ d /i 9 



E 



2fi 



dE n?k dk 



(123) 



and also make use of the relations (which follow from Eq.(9.1.30) of the handbook by M. 
Abramowitz and I. A. Stegun) 



d 

dz 



(124) 
(125) 



where J^x{z) stands for either jx{z) or y\{z). Therefore 



fi d 



K^k dk 



jxjkr) 



A+2 



d 



h?k d{kr) 



j\{kr) 



" mk {krY+^ ^2-^A+ili^,rj 



and thus 



Jx{E,r) 



s^^HEn r) = - (-^ 



dE 

j\+n{kr] 



LIT \ ^ 

-^j jx+n{E,r) 



E=Eo 



1 



fir 



irr 




Ji+n{kr) 



(126) 



(127) 



E=Eo 



As it should be (since j\{E,r) is single-valued), the expansion coefficients si!'\Eo,r) do not 
depend on the choice of the sign of the momentum k^ = ±^/2flEQJW . Indeed, 



3x+n{kr) 



irkr Jx+n+il2{kr) 



nr Jx+n+i/2{kr) 




f^X+n+l Y 2 f^X+n+1 

and since (see Eq.(9.1.35) of M. Abramowitz et al.) 



f^X+n+l/2 



(128) 



(129) 



the numerator and denominator in Eq. (11281) acquire the same phase factor when k changes 
its sign. 



31 



Finding the derivative d^yx{E,r) is a little bit more complicated. The first derivative can be 
written as 



d 



yx{E,r) 



l-A 



^ [ikr)^yxikr)] - -J— [h{k){kr)^jx{kr)\ 



dE ' t9k \ d{kr) '- d{kr 

Using Eq. (11250 and explicit form of the function h{k) given by Eq. ( 1HT1) . we obtain 
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k^-'yx^iikr) - e-'h{k)3x^^{kr) - —k^-'jxikr) 



= -^yx-i{E,r)-—hi{kr) 
where we introduced an auxiliary function 



(130) 



whose derivatives can be found using the following recurrence relation 



(131) 
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I.e. 



dE 
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2nfi ur 

^2 JX,n+l ~r /A— l,n • 



(132) 



Repeatedly using the relations (IISOD and (11320 . we can calculate any number of the derivatives 
d^yx{E,r) needed for finding the expansion coefficients (1122^ . 
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Figure 1: Deformed contour for integrating differential equations (I14II15P and (I22ll23p . 
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-16.2643650096 





-6.2711504590 





0.5036180960 


2 X 10-1^ 


4.9422440057 


0.0000588188 


7.1050168573 


0.5710776714 


7.9987409699 


3.6684977768 


8.5025637363 


7.7605743107 


8.5937554145 


12.3052581635 


8.3193121385 


17.0922638769 


7.6952969586 


21.9663916836 


6.7436612278 


26.9367555304 


5.5244040747 


31.8688591621 


4.0100103640 


36.8118853195 


2.2603614329 


41.6490284540 



Table 1: Spectral points E = Er — iT/2 (in the units 10.96 meV) generated by the 
potential ( !84l ). Their distribution on the complex energy surface is shown in Fig. [3l 
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Figure 2: Domains V for the potential ( !84l) . defined by Eq. (EH) (within the solid curve) and 
by Eq. (!53l) for the rotation angle = O.OStt (within the dashed curve). The energy is given 
in the donor Hartree units (see Sec. E])- 
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Figure 3: Spectral points generated by the potential (I54|) . Their numerical values are given 
in Tabled 
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Figure 4: Model potential (1541) measured in the donor Hartree units 10.96 meV as a function 
of the distance measured in the units of donor Bohr radius 101.89 A. The dashed curve IS a 
typical potential for an empty two-dimensional quantum dot. 
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Figure 5: S-wave scattering phase-shift for the potential ( !84l) . In the curve "A" the sharp 
jumps in tt (corresponding to the first two extremely narrow resonances) are missing and as a 
result it does not obey the Levinson's theorem. 
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Figure 6: Comparison of the exact (thick curve) and approximate (thin curve) S'-wave cross 
section for the potential (1840 . The approximate curves are obtained using the expansion (15^ 
with = 4 (five terms) near the points Eq = 1 x 10.96 meV, Eq = 5 x 10.96 meV, and 
Eq = 7 X 10.96 meV. The insert shows a magnified fragment of the curves near the second 
resonance. 
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Figure 7: The domains within which the Jost function for the potential ( 1840 is reproduced, 
using the first five terms (A^ = 4) of the expansion (169 p . with the accuracy better than 1%, 
5% and 10%. The expansion is done around the point Eq = 7 x [10.96 me V] on the real axis. 
The star shows the third resonance given in Table [H 
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Figure 8: Spectral points corresponding to the third and fourth resonances of the potential 
( l84l) . Stars show their exact locations. Open and filled circles are obtained using the expansion 
( l69l) with N = 2 (three terms) and = 4 (five terms), respectively. The expansion is done 
around the point Eq = (7.55 — 21.06) x 10.96 meV, which is in the middle between these two 
resonances. 
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